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Abstract 

The principles behind the detailed results of a light-front mean field theory 
of finite nuclei are elucidated by deriving the nucleon mode equation using a 
simple general argument, based on the idea that a static source in equal time 
coordinates corresponds to a moving source in light front coordinates. This 
idea also allows us to solve several simple toy model examples: scalar field in 
a box, 1+1 dimensional bag model, three-dimensional harmonic oscillator and 
the Hulthen potential. The latter provide simplified versions of momentum 
distributions and form factors of relevance to experiments. In particular, the 
relativistic correction to the mean square radius of a nucleus is shown to 
be very small. Solving these simple examples suggests another more general 
approach- the use of tilted light front coordinates. The simple examples are 
made even simpler. 
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I. INTRODUCTION 



Light front techniques [|TJ- ||have long been used to analyze high energy experiments 
involving nuclear and nucleon targets. In the parton model ratio k + /p + where k + = k° + k 3 
is the plus-momentum of the struck quark and p + is the plus-momentum of the target turns 
out to be equal to the Bjorken variable XBj- Quark distributions represent the probability 
that a quark has a plus-momentum fraction XBj = k + /p + . In nuclear physics we are often 
concerned with the distribution functions which describe the plus-momentum carried by the 
nucleons within the nucleus. Such distributions, which depend on k + of the struck nucleon, 
are needed to analyze deep inelastic lepton-nucleus scattering, and enter also in the analysis 
of high- momentum transfer nuclear quasi-elastic reactions such as (e,e'), (e,e'p), (p,2p). If 
the light front formalism is used, these distributions are simply related to the square of the 
nuclear ground state wave functions. Many nuclear high momentum transfer are planned, 
so that it is necessary to derive a relativistic formulation which uses k + variables closely 
related to experiments and which incorporates the full knowledge of nuclear dynamics. 

If the relevant nuclear wave functions depend on k + , the canonical spatial variable is 
x~ = x° — x 3 . This leaves x + = x° + x 3 to be used as a time variable, with the light 
front Hamiltonian (x + development operator) as P~ = P° — P 3 . These are the light front 
variables of Dirac. There are clear advantages in using these variables, but a principal 
problem arises because the use of x + as "time" and x~ , x± as "space" involves the loss of 
manifest rotational invariance. This is especially important in nuclear physics because the 
understanding of magic numbers rests on the 2j + 1 degeneracy of single particle orbitals 
of good angular momentum. Our previous papers JlO| showed how light front techniques 



could be used to derive a mean field theory for finite-sized nuclei with results that respect 
rotational invariance. 

There were three features that emerged from that detailed mean field theory treatment: 

• The variational principle which leads to the field equations involves minimizing the 
expectation value of P~ subject to the constraint that the expectation value of P + is 
the same as that of P~ . 

• The meson field equations of the light front (LF) formalism are the same as those for 
the equal time (ET) formulation, except for the replacement z — > —x~/2 in going from 
the ET to the LF formulation. 

• The approximate solution of the LF nucleon mode equation, was a phase factor times 
a solution of the ET nucleon mode equation, again evaluated with z — > —x~/2. 

The purpose of the present paper is to explain these features in a more general way and to 
present some examples of solved problems using light front techniques. The latter is intended 
to build up the reader's confidence in the idea that many different problems can be solved 
using light front techniques. The essential feature that allows a simple explanation of the 
three features is that in each case the mode equation involves external potentials that do not 
depend on x + . We show that such problems are simply related to cases of static potentials 
in the ET formulation. Since using the ET formulation is natural, our basic logic begins by 
noticing that static sources in the usual ET formulation (position (x, y, z) fixed for all time 
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t) correspond to sources moving with constant velocity in the LF formulation (because z + 1 
changes with respect to z — t as t changes. But static sources depend on the variable x + 
in a manner simple enough to be easily removable with a simple transformation to a source 
which is static in light front coordinates. The result is a LF theory in which x + is absent, and 
in which the operator P + is a part of the x + development operator. The transformation is 
applied to the Dirac equation in Sect. II. The same technique is used to solve four problems: 
scalar field in a box, 1+1 dimensional bag model, three dimensional harmonic oscillator and 
the Hulthen wave function in Sect. III. We also consider the resulting light front momentum 
distribution function as a function of p + and p±, and also study the related electromagnetic 
form factors. For the examples we consider, the only difference between our present approach 



and the full calculations in Ref. |10| is in the support properties. In the complete theory 
the wave functions (of positive energy solutions) have no components with k + < 0. This is 
not kept in the present simple approach, and the sizes of errors introduced depend strongly 
on the mass of the particle involved. For nucleons, the errors are shown to be very small 
indeed. 

The transformation to a light front static source suggests the use of a new set of 
coordinates-tilted light front coordinates: "time"r = x° + x 3 , "space" ( = —x 3 . With 
these coordinates a static source in the ET formulation is also a static source in the tilted 
LF formulation, and the requirement of the constrained variational principle that the oper- 
ator P + be part of the "time" development operator is satisfied from the very beginning. 
Tilted coordinates allow us to easily compare LF and ET calculations in problems where a 
single particles moves in an external potential (fixed source, mean field type problems). We 
introduce these coordinates in Sect. IV and demonstrate some utility in Sect. V by applying 
these to solve the simple problems of Sect. III. A brief discussion of the results is presented 
in Sect. VI. 

In this paper, our use of tilted coordinates is restricted to fixed source, mean field prob- 
lems in which one uses quantum mechanics and not field theory. This is done to provide 
the reader with some intuitive ideas about these coordinates and to illustrate their practical 
use. However, this does by no means exhaust the potential applications of tilted coordinates. 
There are many genuine quantum field theory problems that do involve also external fields, 
such as heavy-light systems in QCD (e.g. B-mesons) or a QFT treatment of large nuclei. 
For such systems, the relation between ET and LF is more complicated than the situations 
discussed here because the microscopic degrees of freedom (the quanta) in terms of which the 
theory is formulated are different in these two approaches and the LF approach can provide 
new insights into observables that probe light-like correlations. The use of tilted coordinates 
is that they provide a significant shortcut towards constructing the LF Hamiltonian. 



II. ALTERNATE DERIVATION OF THE NUCLEON MODE EQUATION 

In Ref. flH]] the constrained minimization (of P~) led to a new equation for the single- 
nucleon modes (Eq. (4) in the short paper, (3.28) of the long paper). In order to illuminate 
the physics we present another, more heuristic, derivation of the nucleon single particle wave 
equation. The first step towards formulating the mean field approximation in light front LF 
coordinates is to develop the formalism for static potentials (static in a rest frame!) on 
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the LF. For this purpose, we start from the Dirac equation for a static potential in normal 
coordinates 



i^df, - m - g s V s if) - gvl V {f)} if)' = 0. 



(2.1) 



Since 7 = (7 + + 7 )/2 and since couplings using the "bad" component 7 are difficult to 
handle in the LF framework, we perform the transformation 



where 



<9,r = Vr 



0- 



and T does not depend on time. Using the transformation 
to the vector field only via 7 + 7 s (the LF-7 + component!) 



(2.2) 
(2.3) 

one finds that if)' couples 







'i-fdp + Z7± ■ 0±T) - m - g s V s (?) - g v (l° + 7 3 ) V (f)] $ . 



(2.4) 



Simply rewriting Eq. ( p.l|) in LF coordinates, with 7^ = (7 ± 7 s ) and x ± = (x° ± x 3 ), 
yields 



(ry + <9 + 27 <9 + ) + [dx + g v d±T(f ± 



or — x 



X 



2 

$' = 0. 



— m 



(2.5) 



Even though the potential is static in the equal time formulation, the Dirac equation for 
the same potential in light-front coordinates is LF-"time", i.e. x + , dependent. The physical 
origin for this result is that a static source in a rest-frame corresponds to a uniformly moving 
source on the light-front. Given that the time dependence of the external fields is only due 
to a uniform translation, it should be easy to transform Eq. (|2.5| ) into a form which contains 
only static (with respect to x + ) potentials. For this purpose, we consider the equation of 
motion satisfied by Dirac fields which are obtained by an x + (LF-time) dependent translation 



Using P + = -i2 



dx 



if>'(xj_,x , x + ) = e 
' ' we find 



-P+/2 



if>(x±,x ,x + ) 



(2.6) 



Jx+P+/2 



f 



X 



-IX+P+/2 



f 



—X 



e ix+P+/2Q- e -ix+P+/2 = Q- 

so that the equation of motion for if> takes the form: 
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1 , 1 _ , x~ 

[-27+(<9 - d + ) + -27 d + + «7j_ • (<9j_ + ig v T(f ± , — — )) - m 

-g S V s (r ± , - X —) - gvl + V~(f ± , -^)}iP = 0, 



(2.7) 



(2.8) 



with V~ = Vq. The translated fields satisfy an equation of motion with potentials that 
do not depend on x + . Moreover, the static potentials evaluated at r correspond to light 
front potentials evaluated at (fx, — ^-). A simple derivation of this can be obtained from 
evaluating z = (x + — x~)/2 at x + = 0. 

That the result ( |2.8| ) is the same as the equations for ij)± in Ref. PI . can be seen by 
making a decomposition into a dynamical and a constraint equation. Multiplication of Eq. 
by 7 + from the left yields a constraint equation (no x + -derivative !) 



id + i/j- = ia ± ■ (£>x + igv(dxT)) + (3m + V s ip+ 



(2.9) 



where, as usual, ipx = |7°7 ± V ; - Multiplication of Eq. (|2.8|) by 7 from the left yields an 
equation for 



i(d~ - 8+ - ig v V-)$ + = \iax ■ ± + ig v (dxT)) +(3m + V s ] V>- 



(2.10) 



One may use the constraint equation (|2.9| ) to eliminate ip- in Eq. fl2.10|) to obtain the 



equation of motion for for the dynamical degrees of freedoms. The results (O) and ( |2.10| ) 
are the desired equations. 



III. SIMPLE PROBLEMS 

In order to illustrate the application of the LF formalism for static external potentials, let 
us consider a few simple examples: a scalar field in a box with vanishing boundary conditions, 
a Dirac field in a "bag" with bag boundary conditions, the 3 dimensional (scalar) harmonic 
oscillator, and the 3 dimensional Hulthen potential. 



A. Scalar Field in a Box 

Let us first consider a box of length L (extending from to L) with vanishing boundary 
conditions and determine the eigenstates in the equal time framework. The two degenerate 
solutions to the free Klein-Gordon equation are left and right moving plane waves 4>±(x, t) = 
e -iEt e ±ikz w h ere _£> — yjk 2 + m 2 . Using the boundary conditions to match coefficients and 
wave number in the superposition of these two solutions one obtains the familiar solution 

(f) n (z,t)=e- iEnt sm(k n z), (3.1) 

where 

k n = — - n = 1, 2, 3, ... 

E 2 n = m 2 + k*. (3.2) 

In a light front calculation the boundary conditions depend on x + because z = (x + — 
x~)/2. This dependence can be eliminated using transformation of the form of ( |2.6|) . We 
start with the Klein-Gordon equation 
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(d + d~ + m 2 )4)(x- , x + ) = 0, 



(3.3) 



and make the transformation 



-p+ 



/2 <f> 



to find 



(d + (d- - d + ) + m 2 ) <f>(x~, x + ) = 0, 



with (according to the translation |2.7|) the boundary condition 

0(x" = 0,-2L,x + ) = 0. 

We find solutions of the form 



(x ,x + ) 



e -ik~x+ /2 e -ik+x~ /2 



with the dispersion relation: 



fc = t- + k + . 

k + 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



Before taking the boundary condition into account, there are two linearly independent so- 
lutions e~ %kaX I 2 and e~ lk b x / 2 , each of which has the value of k~ provided by Eq. ( |3.8| ). 
Imposing the boundary condition at x~ = leads to the form: 



(x + , x ) = e 



ik~x+/2 ^ e ik+x~/2 _ p ik+x~/2 



The vanishing of <fr at the other boundary (x = —2L) implies 

n = 1,2,... 



(3.9) 



(3.10) 



which constrains the allowed energy eigenvalues k n . Using Eq. ( |3.8| ) for k a = k b = k~ 
shows that the two independent solutions to Eq. (|3.8|) are related by 



kl 



m 



(3.11) 



and thus the quantization condition for the LF momenta in a stationary box can be written 
as k+ n — {jr- = n^y. Hence the quantized energies satisfy 



m 
k+ 



2tt\ 2 9 
n— J + 4m = 4 



- k + \ + 4m 2 



7T 

n— ) + 



(3.12) 



which is consistent with Eq. (|3.2|) because k n is identified with 2E n . 



The fact that the LF-ansatz Eq. ( |3.9|) consists of two waves — both moving in the same 
direction (both LF-momenta positive!) - has a very intuitive interpretation. In ordinary 
coordinates, the stationary solutions for a particle in a box with hard walls is described by 
a superposition of plane waves moving in opposite directions. Since the LF momentum of a 
plane wave is given by the sum of the particle's momentum in the z-direction and its energy, 
the two waves moving in opposite directions will have different LF-momenta. However, 
both LF momenta are positive, since the ET energy is larger than the absolute value of the 
momentum. 



B. 1+1 Dimensional Bag Model 

As a second example, let us consider the bag model. For simplicity, we will only consider 
the 1+1 dimensional case here. In 1+1 dimensions, Dirac matrices are 2x2 matrices and 
we choose to work in the chiral representation 

»-(!!) (3 ' 13) 

In order to solve the Dirac equation for a "bag" 

i (lodt ~ lid x ) i/j = mi[) (inside) 

n^^ip = iip (at boundary), (3-14) 

where is normal to the surface, we make a stationary wave ansatz 

?P= ( Xl ) e - iEt , (3.15) 



yielding 



X2 



EX2 - 1X2 = m Xi 

Exi + % = m X2 (inside) 

X2 = iXi (x = L) 

X2 = ~iXi (x = 0). (3.16) 



The case m = is particularly easy to solve, and one finds 

Xx = e~ iEnX 

X2 = ~^ lEnX 

En= { U + Dl n = 0,l,... (3.17) 

Turn now to the light front calculation. We use the same representation for the 7-matrices 
as above ( |3.13| ), yielding 

H°o) -r-dl). (3,8, 
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Our version of the Dirac equation is given in 1+1 dimensions (after the transformation (|2.6|) ) 
leads to the 1+1 dimensional version of (12. 81) ) as 



id ip + = id + ip + + rm/>_ 

id + ip- = mip+. (3.19) 
We look for a solution of the form: 

,jj{ x ) = e - tp ~ x+/2 tP(x-), (3.20) 



which gives for ( |3.19| ): 



p ip + = id + ip + + imp- 

id + tp- = mip + , (3.21) 



and therefore 



2 

m 



p>+ = ^— + id + ^ + . (3.22) 

Using plane wave solutions of the form ip + = e~ tk ^ x I 2 or ip + = e~ %k t x / 2 ; leads immediately 
to p~ = jj- + k^ b , so that once again k£ = m 2 /k^. The boundary conditions of Eq. (|3.14|) 

a ,b 

are given in light front coordinates as 

ip- = iip+ {x~~ = —IV) 

ip- = -iip + (aT = 0). (3.23) 
For m = the solutions to Eqs. ( |3.21| ) are given by 

V> + = e~ ip ~ x ~ /2 

= const, (3.24) 

so the boundary conditions Q3.23|) lead immediately to 

Pn = {n+^)— , (3-25) 

which is 2E n . 



2/ L 



C. Three Dimensional Harmonic Oscillator 



In ordinary coordinates, the Klein-Gordon equation for a 3 dimensional relativistic os- 
cillator 

(3.26) 

closely resembles the equation for a nonrelativistic harmonic oscillator and we can thus 
immediately write down its energy eigenvalues 



E 



-Q, , 2 , -*2 

p + m + kx 
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E 2 n = m 2 + u(n + l) , (3.27) 



where uj 2 = 4k. 

The light front version of the Helmholtz equation is 

(d~ d + - V\ + m 2 + k(x\ + (x + - x- f/A)) = 0. (3.28) 

Once again we see a familiar pattern. The static potential in ET coordinates becomes a 
'time" or x + dependent potential in LF coordinates. Once again this dependence is of a 
simple form; it can be transformed away using 



= e- l — X . (3.29) 
The use of Eqns. ( ^.6| ) and ( |2.7| ) then leads to the result 

(V - d + )d + - Vi + m 2 + k(x 2 ± + (^-) 2 )^) X = 0. (3.30) 

In the absence of interactions the P~ operator would consist of the usual term plus an 
additional p + operator. We look for a solution of the form: 

x (x)=e-^ x+ / 2 X (x-,x ± ). (3.31) 
This, along with completing the square, leads to the result 

(~(d + + ^P~/2) 2 - Vi + m 2 + k{x\ + (^) 2 ) X = fc/2)V (3.32) 

One converts the operator (d + + ip~/2) 2 to (<9 + ) 2 using yet another transformation: 

X(x~,x±) = e- ip - x ~ /4 F(x-,x ± ) (3.33) 

to find 

(-(d+) 2 - Vi + m 2 + k(x 2 ± + F = (p-/2) 2 F. (3.34) 

This is the same form as the equation in the equal time coordinates and p~/2 takes on the 
values of E n . 

We seem to be getting the same results as in the ET development. So one might wonder 
why we are doing the light front at all. The point is that we are able to compute the light 
front wave functions that depend on x~ , or in momentum space depend on p + . The wave 
functions of the ground state is given by 



1 x~ 2 



Xo 



(x~,x ± ) = e ip o*-/ 4 jVoexp + ~r) ■ (3.35) 
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The number density no(p±,p + ) is defined as the square of the momentum space version 
of Xo- This quantity is accessible in high energy proton and electron nuclear quasi-elastic 
reactions. It is useful to define the light front variable 

«^ + /(f)- (3-36) 
Then one may easily determine that 



no(p±,a) = N e ^ . (3.37) 

Note that one finds the same p± distribution for each value of the variable a. This is not a 
general feature of light-front wave functions, as we show in the next sub-section. 

In an exact calculation, the number density should vanish for values of a that are not 
between and 1. This is referred to as the support, which we examine here. How large can 
the value of a be? For large values of the particles mass m, (m ^> k 1 / 4 ) (which corresponds 
to the situation of nuclear physics in which the product of the nucleon mass and the nuclear 
radius is a very large number) the value of a must always be close to unity. If m — > then the 
behavior is ~ e~ 3(a_1 ) . We may better understand this factor by considering an extremely 
simple model- the nucleon consists of three massless quarks moving in an harmonic oscillator 
potential. We shall compute the quark distribution function q(x), which here is the same as 
n integrated over p± and evaluated as a function of the Bjorken variable, x. For massless 
quarks, Eq. fl3.27|) gives 



( y) 2 = 3v^- (3-38) 

Furthermore, in the target rest frame P£ is the mass of the nucleon 3-£-. Thus using 
Eq. G3.36Q a = 3p + /P£ = 3x. The last equality is from the parton model in which the ratio 
of the quark and target plus momenta is 2 m N v = x ~ Thus, we find 

q(x) oc e- 27 ^- 1 ^ 2 . (3.39) 

In this limit of massless quarks, the value of x (in this un-evolved) quark distribution function 
is constrained to be very close to 1/3. One would naively expect the value of x to easily 
exceed unity, since we have used a mean field model for a three-quark system. This does 
not occur. If x — 1 the value of the q is e -12 its peak value, which is reasonably small. 
If x approaches 0, one finds a factor e~ 3 , instead of the required so there is an inherent 
inaccuracy of some 5%. If one included a non-zero value of the quark mass, the support 
properties would be improved because 

g(x)oce- ( * +27)(a; - 1/3)2 . (3.40) 

In constituent quark models m? = yfn, so that at x = one finds a factor of e~ 4 . In nuclear 
physics ^= is of the order of (5-R^/Fm) 2 which taking = 4 Fm, yields e -400 at x = 0. In 
that case, there is no problem with the support. 
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It is of interest to compute the electromagnetic form factor of the ground state. This 
quantity has been often measured in elastic electron scattering, and the sizes of nuclei have 
been determined as one of the classic achievements of nuclear physics. Interest in this topic 
has been revised because of a recent proposal to Jefferson Laboratory to use parity- violating 



electron scattering to measure the neutron radius JTT|. A high precision is needed and can 
be obtained provided one knows the proton distribution. Therefore one needs to examine 
the influence of small effects such as relativistic corrections. One works using a reference 
frame in which the plus component of the four vector q^ of the virtual photon vanishes, 
so that Q 2 = —q 2 = q]_. In this case the form factor F(Q 2 ) (matrix element of the plus 
component of the electromagnetic current operator) is given by 

F{Q 2 ) = J d 2 p ± da Xo(a,P±)Xo(a,P± - aq±), (3.41) 

in which the influence of relativity appears in the integral over a and the factor a. For the 
harmonic oscillator ground state we find: 



.2 



F(Q 2 ) = N J d 2 p±daeV K ^ J e'^^ , (3.42) 

where 

_2 

^- = m 2 + 3v^. (3.43) 

Our purpose here is the study of nuclear physics, so we are interested in the non- 
relativistic limit and the corrections to it. To this end, we define a variable p z using 

a = 1 + — . (3.44) 
m 

_2 

The non-relativistic limit of ( |3.42| ) is obtained by letting m approach infinity. Then = m 2 
and we find 



F NR (Q 2 ) = N NR [ 

J m 



e (3.45) 



The mean square radius — & dF }q2 - |q 2 =o> is given by 

R% R = 3 -±. (3.46) 

The leading corrections to this will be of order p 2 /m 2 ~ ^Jk/wl 2 . We define a semi- 
relativistic limit SR via the use of (|3.42|) and keeping the leading correction terms. Per- 
forming the straightforward evaluations leads to the result: 

6 = Wr = (3 ' 47) 

or 



This corresponds to very small (0.004) effects for large nuclei A ~ 200. 
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D. Light Front Hulthen Wave Function 



Any static potential of the form V(x 2 = x\ + z 2 ) can be solved on the light front. The 
transformation ( |2.6| ) corresponds to including the p + term in the x + development operator 
and a simple prescription of replacing z by —x~/2 in V. We present here the solution for 
the Hulthen potential. This allows us to demonstrate an interesting contrast between the 
implications of different forms of potentials. In the equal time formulation we have the wave 
equation: 



E 1 



f + m 2 + V H (x 2 ) 



in which V H is chosen so that the lowest energy solution is 

0(r) = N(e~ ar - e- br ), 
where b > a. The eigenenergy is given by 



E = V: 



— er 



The light front version of the wave equation is 



Pn 



\\X , X± / 



-21 



d 



+ P 'i + v H [x 2 ± + 



dx~ 



x 



X(x ,x±). 



(3.49) 



(3.50) 



(3.51) 



(3.52) 



The lowest value of p n /2 is clearly the same as E of Eq. ( |3.51| ), and the wave function is 
given by 



,fi) = e ip " x ~ /4 N^ 



exp —a\ x ± + 



x 



exp — b\ x\ + 



X' 



,_2 



(3.53) 



The momentum distribution \p±,p + ) obtained here provides an interesting contrast with 
that of the harmonic oscillator (|3.37| ). We find 



N, 



H 



(a 2 +p 2 ± + (p-/2) 2 ( a -r 



((b 2 +p 2 ± + (p-/2) 2 (a-l) 



(3.54) 



It is clear that one finds a different p± distribution for each value of a. The distribution 
is a broader function of p± for larger values of a; see the Fig. |l[ The results in the figure 
are obtained using m=.94 GeV, E=.932 GeV and b = 5a. An experimental hint of such a 
behavior has been found recently [p^] . 



12 



FIGURES 
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FIG. 1. Light front momentum distribution as a function of a and p±. 



IV. TILTED LIGHT-FRONT COORDINATES 



We have worked out several examples involving potentials that are independent of time 
in equal time coordinates. This leads to an x + -dependent interaction when light front 
coordinates are used. However this dependence is very simple, and in each case we have 
removed this by using a version of Eqs. (|2.6| ). In each case we have derived a light front 
wave equation in which the kinetic energy has the same form as the standard LF-kinetic 
energy plus an additional term linear in the momentum which is identical to the recoil term 
in the static source formalism on the LF! It is worthwhile to see if there is a more general 
way to remove this dependence, once and for all, by finding a set of coordinates in which 
static sources in equal time coordinates are also described by static sources in a coordinate 
system that is very much like that of the light front. 

For this purpose, we introduce new "tilted" LF-coordinates 




(4.1) 



i.e. 



X° = T + ( 

x 3 = -C 



(4.2) 



with x± as usual. These coordinates very much resemble LF-coordinates, since 
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FIG. 2. Comparison between lines of equal x° and x s (full lines) and lines of equal r and C 
(dashed). 

• surfaces of constant r are the usual LF hypersurfaces and thus quantization is very 
similar to LF quantization. 

• correlation functions at r = in the ( direction yield the usual LF distributions. 
Therefore, a lot of the familiar LF-phenomenology (e.g. structure functions and wave 
functions) can still be used. 

However, in contradistinction to light-front coordinates, static sources in a rest-frame are 
also described by static sources in above tilted coordinates, since 



d = d T 

d 3 = d T - <9 C , 



(4.3) 



i.e. if OqV = then d T V = 0. For a static (dtV = 0) potential, one thus finds for the static 
potential V(x±,x 3 ) 



V(x ± ,x 3 ) = V(x ± ,-C) 



(4.4) 

Furthermore, using Eq. ( |4.3| ), one finds for the longitudinal part of the kinetic energy 
operator 

dl - dl = 2d T d c - d 2 c (4.5) 
The Lagrangian for a scalar field interacting with an external field can then be written as 

1 r. 



C 



d, 



- (m 2 + V{x)\ (f) 



(4.6) 
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2 m 2 + V{x) 



in tilted coordinates, the canonical momenta are the longitudinal "space" -derivatives of the 
fields 



dC 

^w^r^^ 9 -* (47) 

— just like in genuine LF-coordinates. This should not come as a surprise, since the constant 
"time" surfaces in our tilted coordinates and in LF coordinates are identical (Fig. |2|). 

The kinetic energy operator differs from the one in LF-coordinates. Using Eq. ( |4.5|) in 

the form 2£V — V 2 — P ± — m 2 = 0, where £ is conjugate to r and V is conjugate to (, one 
finds the kinetic energy of a free particle in tilted coordinates: 

£ = ^ + P. (4.8) 
For non-interacting systems, this is identical to the expression derived in LF-coordinates 



including the "Lagrange multiplier" [ID]. 



The point of using the tilted LF coordinates is that the term V is included automatically. 
The energies £ obtained using tilted LF coordinates are the same as P~ obtained in the 
usual LF dynamics including the Lagrange multiplier term. We expect that the £ of tilted 
LF must be the same as the P~ of LF because both formulations use the same quantization 
(i.e. equal 'time') hypersurfaces. The only place where they differ is the space direction, 
where the difference between the spatial coordinates x~ and £ is a shift x°. This is a simple 
translation — hence the term linear in the momentum in the tilted coordinates Hamiltonian. 
However, it turns out that in the mean field approach one needs to introduce a Lagrange 
multiplier term to fix the momentum, and the two Hamiltonians have an identical form. 

We shall see how this works out by solving the same simple problems as above using the 
tilted LF dynamics. 



V. SIMPLE PROBLEMS ON THE TILTED FRONT 

In order to illustrate the application of the LF formalism for static external potentials, 
let us consider a few simple examples: a scalar field in a box with vanishing boundary 
conditions, the 3 dimensional (scalar) harmonic oscillator and a Dirac field in a "bag" with 
bag boundary conditions. 



A. Scalar field in a box 



In order to derive the tilted LF solution for this example, we first note that vanishing 
boundary conditions are frame independent, i.e. a stationary box with vanishing boundary 
condition in an equal time framework corresponds also in tilted LF coordinates to a sta- 
tionary box with vanishing boundary condition. Before taking the boundary condition into 
account, there are two linearly independent solutions e~ tka< * and e~ lk <> for a given energy, 
which is the same according to Eqs. (|4.8| ) as that of Eq. Q3.SQ. 
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The only superposition which satisfies the requirement that the wave function vanishes 
at ( = is given by 

= e - ik « c - e~ ik ^. (5.1) 

Vanishing of <f)(x~) at the other boundary = —L) implies 

2n 

K-kt = n— n = l,2,... (5.2) 



which constrains the allowed energy eigenvalues k n . Using Eq. (|4.8| ) for k a = k b = k r 
shows that the two independent solutions to Eq. (|3.8|) are related by 



2 

m 



K = ^ (5-3) 

and thus the quantization condition for the LF momenta in a stationary box can be written 
as kan ~ = n T"- Hence the quantized energies satisfy Eq. (|3.12j ), which is consistent 
with Eq. ( |3.2|) provided one identifies k~ with 2E n ( |4.3|) . 



B. 1+1 Dimensional Bag Model 



Making a stationary wave ansatz 

4 = 



1>- 



-i£r/2 



(5.4) 



one thus writes the equation of motion for a free Dirac particle ( [3.14 ), in tilted coordinates 
by multiplication by either 7 + or 7". This gives 



—id^ip- 



mi/j- 



(5.5) 



where = 4^. We use the second of Eqs. ( |5.5|) in the first to to obtain the "dynamical" 



7 C - d( 
component: 



m , -a 
10 t 



(5.6) 



As is the case for free scalars, one finds in general two linearly independent plane wave 
solutions for each energy ip + = e~ lka< * and ip + = e _ifc f> where £ = pp- + k^ b and k b = —p. 

a,b ' ' k °- 

Note that the boundary condition ( 3.14Q mixes ip + and ip- 



(C = -L) 
(C = 0). 



(5.7) 



This mixing should not come as a surprise, since the boundary condition arises from assuming 
an infinite mass for the fermion outside the bag. 
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For m = the solutions to Eqs. (|5.5| ) read 



and thus from the boundary conditions 

£ = 



const. 



1\ 2tt 



n 



(5.8) 



(5.9) 



Using again S n = 2E n , we find that the spectrum obtained in tilted LF coordinates is 
consistent with Eq. (|3.17|) . 



C. Three Dimensional Harmonic Oscillator 



In order to solve the Helmholtz equation in tilted coordinates 



V£4> = [m 2 + k (x 1 + C 2 ) + f± + V' 



we first complete the square, yielding 
£ 2 ( £\ 2 



m 



(5.10) 



(5.11) 



First of all, we note that a shift in V, included by multiplying the wave function by a 
suitably chosen phase factor can be made to absorb the dependence of the r.h.s. of Eq. 
(|5.11 ) on S. Furthermore, the transverse and longitudinal dynamics in Eq. (|5.11| ) completely 
separates, yielding 



^ = m 2 + uo L (n L + ^ + lu ± (n ± + 1) 



(5.12) 



where uul — — 4k = lo. Taking into account that S n = 2E n we thus find that the 
relativistic harmonic oscillator formulated in an equal time framework and in tilted LF 
coordinates yield identical spectra. 



D. Dirac Equation with Static Potential 



We express Eq. (|2.5| ) in tilted coordinates to find: 











V s (f±, -C) - J + V (?x, -C)l V' = 0. 



rn 



(5.13) 



The potentials V$ and Vq are functions of r\ + ( 2 , so we expect rotational invariance to be 
immediate. Indeed, the above is just the standard Dirac equation, except for the appearance 
of 7 + instead of 7 . This can immediately be removed by using the transformation of 



Eq. (3.44) of Ref. |T(J. Thus the spectra of the Dirac equation in tilted coordinates is going 
to be the same as that of the usual equal time formulation. 
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VI. SUMMARY 



The three features discussed in the introduction emerge naturally from a more general 
notion that static sources in the usual ET formulation become static sources in the LF 
formulation if one makes an appropriate transformation (|2.6|) . The resulting LF formulation 
for scalar interactions can be summarized as simply including a p + in the x + development 
operator, and evaluating the potentials using the replacement z — > —x~/2. The resulting 
simplicity allows us solve any problem involving a spherically symmetric static source. An 
application of this in Sect. V.D enables us to obtain a momentum distribution which increases 
it width as a function of p± for increasing values of p + . A hint of such behavior has recently 



been observed 12 



We note that the spectrum condition-the feature that light front mode functions have 
support only for positive values of plus-momentum- has not been maintained in any of the 
solutions we present here. The desire to maintain this consistency with the original field 
theory led us to implement a new numerical procedure in Ref. flit)] . In the case of the 
harmonic oscillator, we see that for the general conditions of nuclear physics, the spectrum 
condition is maintained, even though the solution procedure ignores this condition. The 
detailed study of the how important maintaining the spectrum condition, for models other 
than the harmonic oscillator, is will be a topic of future investigation. The numerical results 
of Ref. ( [fLOU), for nucleons in nuclei, strongly indicate that there is no strong need to 
maintain the spectrum condition. The presence of a 5% effect for q(x = 0) indicates that 
the maintaining the spectrum condition is important for massless quarks. 

The simplicity of the results for the simple modeled considered, along with the need to 
automate the transformation procedures, suggests the use of a new set of coordinates- tilted 
coordinates in which the p + term is in from the beginning and in which static sources in 
the ET formulation are also static sources in the tilted LF formulation. For the case of the 
simple models considered here, this new formulation is actually easier to apply than the 
usual LF formulation. Whether or not this simplicity survives more detailed problems is a 
matter for future investigation. 
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